Crossover scaling from classical to nonclassical critical behavior. 
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We study the crossover between classical and nonclassical critical behaviors. The critical crossover 
limit is driven by the Ginzburg number G. The corresponding scaling functions are universal with 
respect to any possible microscopic mechanism which can vary G, such as changing the range or 
the strength of the interactions. The critical crossover describes the unique flow from the unstable 
Gaussian to the stable nonclassical fixed point. The scaling functions are related to the continuum 
renormalization-group functions. We show these features explicitly in the large-iV limit of the 
0(N) (j) 4 model. We also show that the effective susceptibility exponent is nonmonotonic in the 
low-temperature phase of the three-dimensional Ising model. 
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Motivated by various experimental results (see e.g. 
Refs. there has been recently a revived interest in 

understanding crossover phenomena driven by the effec- 
tive range of the interactions. If the interactions have a 
finite range R, in the limit in which the reduced tempera- 
ture t goes to zero, the system shows the standard short- 
range nonclassical behavior. According to the Ginzburg 
criterion § this occurs when !«G where t and G are 
respectively the reduced temperature and the Ginzburg 
number. On the other hand, in the opposite limit t ^> G 
the system shows a classical Gaussian behavior. In the in- 
termediate region one observes a crossover between these 
two behaviors. From the point of view of the Wilson 
renormalization-group (RG) theory, this crossover phe- 
nomenon is generally explained by the competition of two 
fixed points: the Gaussian fixed point and the nonclassi- 
cal fixed point that determines the asymptotic behavior 
in the neighborhood of criticality. 

These crossover phenomena are of great importance 
for the understanding of critical phenomena occurring in 
physical systems (see e.g. Refs. ) . Fisher [0 dis- 

cussed experiments on micellar solutions [|l] and argued 
that the apparently nonuniversal results of the critical 
exponents may be explained in terms of a crossover be- 
havior driven by the effective range of the interactions. 
Crossover phenomena are also observed in experimen- 
tal data for the susceptibility of fluids and liquid mix- 
tures ||, and in polymer melts ||. Some understanding 
of the crossover problem is provided by field-theoretic 
calculations (see e.g. Refs. |7|-[l5|). 

The most important issue concerning crossover phe- 
nomena is whether one can define scaling functions that 
are universal and that describe the crossover between the 
classical (Gaussian) and the nonclassical behaviour. In 
order to give an answer to this question one must clarify 
the kind of crossover one is considering. Varying t, while 
keeping the Ginzburg number G fixed, one observes a 
crossover between the critical and the noncritical behav- 
ior. This is obviously non-universal and for t < G it is 
described by the nonuniversal Wegner expansion ]l6| . As 
pointed out by Bagnuls and Bervillier JlTj (see also |18|) 



a universal behaviour can only be obtained in a prop- 
erly defined critical limit. In other words, a universal 
behaviour can be observed only if t <C 1, i.e. the cor- 
relation length is large, in the whole crossover region. 
The problem was properly formulated by Luijtcn, Blotc 
and Binder Jl^,0 , who argued that a universal crossover 
description can be obtained if one considers the simulta- 
neous limits t — ► 0, R — > oo keeping the product be- 
tween t and an appropriate power of R fixed. A Wilson 
RG analysis indicates that this limit is non trivial 
and interpolates between the mean-field and the standard 
short-range behaviour. These ideas have been confirmed 
numerically in the two-dimensional Ising model j2Cj . We 
will refer to the above limit as critical crossover. Thus the 
critical crossover is a crossover from the critical nonclas- 
sical behavior to the critical classical behavior. Extend- 
ing the arguments of Ref. [ p0[ , we define the appropri- 
ate crossover limit in the whole (t, h) plane introducing 
a magnetic Ginzburg number Gh, such that the system 
shows classical behaviour for h ^> Gh and the standard 
short-range behaviour in the opposite case. 

The main point is that in the critical crossover region 
the key role is played by the Ginzburg number G: the 
range of the interaction is only one of the possible micro- 
scopic parameters controlling G. Any other mechanism 
leading to a change of G can give rise to the same critical 
crossover when the appropriate limit is considered. We 
show that the critical crossover driven by R can be re- 
produced starting from a standard (f> 4 theory with short- 
range interactions and taking an appropriate limit of the 
theory when the bare four-point coupling goes to zero. 
The critical crossover functions expressed in terms of the 
renormalized coupling are related to the standard con- 
tinuum RG functions. 

The several open questions on the crossover behavior 
call for a theoretical laboratory where the various con- 
jectures can be verified analytically. The O(N) vector 
model in the large- N limit is ideal for this purpose. In- 
deed, although it maintains many nontrivial features of 
the theory, it allows us to perform exact calculations, 
and therefore an exact verification of the conjectures 
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on the crossover phenomena. Starting from a lattice 
O(N) model with long-range interactions, we calculate 
the crossover functions for 2 < d < 4, and discuss their 
universality in the critical crossover region. We derive the 
equation of state that provides a complete description of 
the crossover region. 

For the sake of definiteness, we consider the d- 
dimcnsional lattice model defined by the Hamiltonian 

U = £ \J(x { - Xj) [<j){xi) - fiixj)} 2 

i 

where 4>{xi) are ^-dimensional vectors. The spin-spin 
interaction J(x) has finite range R defined by 



r2= 1 J2x X J (x) 



2d £.J(£) 



(2) 



The specific form of J(x) is irrelevant for our discussion. 
The normalization of J(x) is chosen so that its Fourier 
transform II = H(k, R) has the low-momentum behav- 
ior k 2 + 0(k 4 ). The Ginzburg criterion || applied to 
the model (|l|) tells us that the theory has a nonclassical 
critical behavior when 



t 



r c «Goc U 2 /( 4 " d ). 



(3) 



In order to study the long-range limit, it is convenient 
to perform a field rescaling with a corresponding rescaling 
of the Hamiltonian parameters: 



R- 



t = R t 



2+ T7T - 



R%, h = Rh. 



(4) 



Keeping R finite, the critical behavior for I — > and 
h — > is nonclassical. For example in the large- N limit 
the critical exponents are 7 = 2v = 2/{d — 2), rj = 0, 
= 1/2, S = (d + 2)/{d - 2). When R -► 00 one ex- 
pects a mean-field critical behavior, characterized by the 
exponents 7 = 2u = 1, 77 = 0, f3 = 1/2, S = 3. This 
change implies a singular dependence on R of the criti- 
cal amplitudes, which has been derived in Refs. pl| , |l9| . 
For example the asymptotic behaviors of the magnetiza- 
tion, of the magnetic susceptibility and of the correlation 
length, for t — > (with h = 0), are expected to be [pl|Jl9| 



(5) 



f = ~2d^R? E X *®° ■ ^ K r*'^- 8 -)/^, (6) 
M=(fi<xfB*W- 1 M*-' t ). (7) 

Moreover, usingthe Wilson renormalization-group ap- 
proach of Ref. 19 1, we find that at t = 



Moct /S R d <- 3 / 5 - 1 ^- d \ 



(8) 



In order to describe the critical crossover from the non- 
classical to the classical behavior as driven by the range 
of the interaction, i.e. keeping u fixed, it is convenient to 
introduce a new Ginzburg number associated with t p9[ : 



G = R 



-2d/(4-d) 



(Nu 



,2/(4-d) 



tx R 2 G 



(9) 



(we have introduced the iV-dependence only to make the 
large- TV limit more transparent). Thus the comparison 
of t with G oc U _2rf /( 4_d ) tells us whether the critical 
behavior is nonclassical (t <§; G) or classical (t » G). 

Following Refs. |l9|,^0|], one may introduce the rescaled 
reduced temperature 



t = i/Goct/G, 



(10) 



and consider the limit t — > and R — » 00 keeping t fixed. 
When i — ► (i — ► 00) the nonclassical (classical) critical 
behavior should be recovered. Extending the RG analysis 
of Ref. {n| to the line t = 0, we also introduce a magnetic 
Ginzburg number Gh oc u( d + 2 )/[ 2 ( 4 ~ <1 )] and 

G h EE R-™/^-*) (JV¥) (d+2)/[2(4-d)] ^ i j Gfc< (n) 

Gft tells us, in the presence of a magnetic field, in which 
regime we are: nonclassical when h <§C Gh and classical 
when h 3> Gh- Correspondingly we define a rescaled 
magnetic field 



h I cx h /G h , 



(12) 



and study the behavior of the theory when h — > 0, i? — > 
00 with /i fixed. The nonclassical (classical) behavior is 
obtained in the limit h — ► (/i — * 00). 

The relations (||) and (|J) suggest the following scaling 
behaviors in the critical crossover limit for h = 



X = xG«iyi)cx X G, 
?=f G«F C 2(F)oc£ 2 G. 



(13) 
(14) 



From Eq. (H) one obtains 



M = M G / Gh ~ F M (h ) oc MG/G h for i = 0. (15) 

F x (x), F^(x) and Fm{x) are expected to behave as 
F x (x) ~ a;" 7 , F^(x) - x~ 2l/ and Fm(x) ~ a; 1 /- 5 for 
x — > 0, and -F x (a;) ~ a; -1 , -^£2(0;) ~ a; -1 and Fm(x) ~ 
2.1/3 f or x _^ OQ _ rp^g corrections to these asymptotic be- 
haviours should be controlled by the corresponding lead- 
ing correction-to-scaling exponents A (see e.g. Ref. 0] 
and references therein). 

It is crucial to notice that in the crossover region the 
relevant new scale is provided by the Ginzburg number 
G, and the critical crossover limit can be expressed in 
terms of G (or Gh) only, i.e. without the explicit use of 
R. The range of the interactions represents a physical 
way to vary G according to Eq. <M), but it is not the only 
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way. The critical crossover scaling functions are therefore 
expected to be universal with respect to the microscopic 
ways one uses to control and vary G. In other words, 
the critical crossover describes the unique flow from the 
unstable Gaussian to the stable nonclassical fixed point. 
Starting from a </> 4 short-range theory, for instance from 
the Hamiltonian (|l|) with R = 1 fixed, one may use the 
bare four-point coupling u to vary G according to Eq. (|^). 
One then recovers the critical crossover behavior in the 
limit u -> 0, t -> with t = t/G cx tu~ 2/{A ~ d) fixed. In 
the context of the statistical approach to polymers the 
limit u — > keeping t cx tu^ 2 ^ 4 ^^ fixed is essentially 
equivalent to the so-called two-parameter model (see e.g. 
Refs. |j22"|-^4"[ and references therein). In this limit the 
crossover functions can be computed in the standard con- 
tinuum </> 4 theory (| |ll],[l^] . A dimensional analysis shows 
that (using the subtracted bare mass and removing the 
cutoff) finite results can be obtained in terms of the di- 
mensionless variable u/t 2 ~ d / 2 = t d / 2 ~ 2 , and no further 
limiting procedure is required. It is useful to introduce 
a renormalized coupling g. Changing variables from t to 
g, one may then show that the critical crossover func- 
tions expressed in terms of g are related to the standard 
continuum RG functions. They describe the physics of 
strongly correlated systems in the whole range between 
the classical and nonclassical critical point, in terms of 
a single physical parameter measuring the ratio between 
the interaction scale and the correlation scale. The criti- 
cal crossover functions for physically interesting systems 
are well studied |) 11 1^] in the fixed-dimension expan- 
sion when d = 3. 

To check explicitly these ideas, let us consider the 
large- TV limit. Before any rescaling the following saddle- 
point equations hold 



„ r 9 6 . „ T f d d k m 2 

M + -(t - to 2 ) = N / - — -j— ^, 

u ' J (27r) d n(n + TO 2 )' 

h 

M= — , 



(16) 
(17) 



where M = (</>) is the magnetization and £ = 1/m a 
dynamically determined length scale. For h = 



X = N/i 



t2 ™2 _ D 2 2 



to = Km. 



(18) 



The critical crossover function F x (t ) can be obtained by 
rewriting Eq. ( |l6| ) for h = in terms of t and x, and 
taking the limit R — > oo (thus G — > 0) with u fixed: 



~ Nu 

t = Nx' 1 + — lim 

O R— >oo 



d d k 



Nx 



(2ir) d n(n + GNx~ 



(19) 



where II = R 2 H. An analysis of the integral in Eq. jl9| ) 
shows that, for 2 < d < 4, the limit exists. It depends 
only on the following property of II: 



lim 



lim IL(q/R,R) 

R^OQ 



(20) 



where boo is a non- vanishing constant. Its explicit value 
is not relevant, since a change of this constant can be 
reabsorbed in a change of normalization for t. Therefore 
our results are universal, apart from a rescaling of t and 
X, for a large class of Hamiltonians satisfying the above 
condition. From Eq. (|l9|) one finally obtains 



-Ml 

N 



K = 1 



K + La 



Nu 



N 



2-d/2 



(21) 



6&L 



F(l - d/2) 
6(47r) d / 2 ■ 



Eq. ( ^l| ) is well defined in the large- A'' limit after proper 
rescalings in N of the fields (i.e. cf> — > yN 4>) and cou- 
plings (i.e. u — > u/N). The expected asymptotic behav- 
iors are reproduced: 



F x (t) (t- 2 /( d - 2 \l + Cl t A - + ...) 



N 



t- 1 {l + c 2 t-^ ■ + ...) 



ioT l 
for t 



0. 

oo, 



(22) 



where A sr = (4 - d)/(d - 2) and A g = (4 - d)/2 are 
the leading correction-to-scaling exponents related to the 
nonclassical short-ranged and the Gaussian fixed point 
respectively. The coefficient of the t A " T correction in the 
expansion around t = (and in general also the other co- 
efficients of the expansion) can be obtained by performing 
the appropriate limit of the nonuniversal Wegner expan- 
sion. We stress that the dependence on the bare cou- 
pling u in Eq. ( f2l| ) can be eliminated by a rescaling and 
a redefinition of t. Moreover, the same equation (modulo 
the above mentioned rescalings) can be obtained starting 
from the A^-vector (nonlinear a) model. 

The analysis of experimental data in the crossover re- 
gion is usually performed by introducing effective critical 
exponents. One can define 7 c s by the logarithmic deriva- 
tive of F x (t). In three dimensions 



7eff 



d\nF x 
dint 



i + (i + c 7 T) 



-1/2 



(23) 



where 



AK/L 2 . 7off(i) is universal apart from a 



trivial rescaling of t. Analogously one may define v e f[ 
and find 



feff 



ldhxF C 2 



Teg 

2 



(24) 



2 dint 

From Eqs. (16) and (pj) one can derive an equation of 



state relating the rescaled variables t, h, and M in the 
critical crossover limit. Simple calculations, involving the 
same integral of Eq. (091) , lead to the equation 



M 2 ~ T h 

h< = 

6N M 



V Af . 



d/2-l 



(25) 



which turns out to_be universal apart from trivial rescal- 
ings of t, h, and M. Moreover it reproduces the correct 
asymptotic behaviours: 
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h~ M s (l + ciM- 1//3 



1 + 



(26) 



where in the nonclassical limit (t — ► and h — > 0) 5 = 
(d+2)/(d-2), /? = 1/2, and p = 2/(d-2). In the classical 
limit 5 = 3, = 1/2 and p = 1. Setting t = in Eq. (H|, 
one can derive the crossover function FM{h) defined in 
Eq. (14). In three dimensions the corresponding effective 
exponent 5 e g is given by 



<iln h 
d\nF M 



2 1 



Oy F M 

6 N 



-1/2 



(27) 



We can also consider the short-range version of the 
model and show that the function F^(t ) = xiNu) 2 ^ 4 ^^ 
satisfies Eq. (|2l]) with boo — > oo. The same arguments 
apply to all other critical crossover functions. This con- 
firms that the critical crossover functions are universal, 
i.e. independent of the mechanism driving G. 

In the symmetric phase, we can define the zero- 
momentum four-point coupling g as 



37V 



9 



Xa 



N + 2 X 2 S, 



led ■ 



(28) 



where xa ' IS the connected four-point correlation function 
at zero momentum. In the large- N limit |2q] 



Ng 



where 



■d-A 



d d k 



1 



Nu 

<i ./ (2n) d (TT+ Gm 2 ) 2 



(29) 



£ 2 . In the critical crossover limit the inte- 



gral depends only on boo- We obtain 



9(0) 



1 



1 



c m 



4-d - 



where Ng(0) = 6(An) d / 2 /T(2 - d/2) = Ng* is the non- 
classical critical value of g, and c g — 2KJ (d — 2)Ld- The 
effective critical exponent associated with g(fh) is 



dhig(m) 
rflnrn 



(d-A) 



1 



9_ 

9* 



(31) 



In the same limit Eqs. fl23j), (124) and (30) imply also 



7eff(ff) = 2v cS (g) = 1 + 



4-d g 
d-2g* 



(32) 



Eqs. ( |3l| ) and ( |32] ) are now independent of u and boo- No- 
tice that only in the nonclassical critical limit (m — ► 0), 
V J cff (9) ~> and therefore the corresponding hyperscaling 
relation is satisfied. This fact is not unexpected because 
hyperscaling is not satisfied at the Gaussian fixed point, 
where g ~ (T - T c )( 4 ~ d )/ 2 for T -► T c . 

One can now easily verify that ipos(g) — (3(g)/ g where 
(3(g) is the Callan-Symanzik /3-function in the continuum 
4> A theory. Analogously the high-temperature exponents 



^ 1.1 




FIG. 1. Effective susceptibility exponent as a function of t 
for the high- (70ft) and low- (7^) temperature phase of the 
three-dimensional Ising model. Here t is the reduced temper- 
ature for the model with Hamiltonian ([^j]) and R is defined 
in Eq. (§). 



lcs(g) and v e s(g) are related to the standard RG func- 
tions 7(g) and v(g) (see e.g. Ref. jy]]) through the rela- 
tions Eq] 



Jeff(g) _ 7(5) 



Vett(g) 



dg 



i(g) - 7ofi(ff)- 



(3°) In the large- TV limit 



2 1 - 1 (g) = v~ 1 (g) = 2 + {d~A)g/g* 



(33) 
(34) 

(35) 



Using the field theoretical approach p ]jHJl2|| , one can 
compute the effective exponents in three-dimensional 
O(N) models. Results for the high-temperature phase 
are reported in Refs. ^JlTf. We extended the compu- 
tation |27j| to the low-temperature phase computing 7^ 
for N = 1. The resulting curves are plotted in Fig. [I]. 
We stress that, apart from the small error of the resum- 
mation procedure — it should be well below 1% — the 
curves in Fig. ^ represent the universal critical crossover 
exponents. Thus experimental and numerical data in the 
crossover region should approach these curves in the ap- 
propriate limit (modulo a rescaling of t ) . In this perspec- 
tive it is possible to understand the lack of universality 
of the results of Ref. || : universality is recovered only in 
the limit u — > 0. For finite values of u one expects cor- 
rections to scaling that eventually disappear as u — > 0. 
Indeed the comparison with the experimental data for 
fluids and liquid mixtures || improves as the effective 
parameter u decreases. 
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Measurements in the critical crossover limit are dif- 
ficult and accurate results are scarce. Experiments on 
micellar solutions |jj observed exponents that were very 
far from the expected Ising values. The exponent 7 was 
even lower than the classical value 7=1. Fisher Q 
interpreted the data as a crossover effect, suggesting a 
standard scaling description. In order to explain the 
data, this interpretation would require j c g to be non- 
monotonic in the symmetric phase. But, as Fig. |l| shows, 
this is not the case for the critical crossover function j e g. 
Therefore, as already observed by Bagnuls and Bervil- 
lier |50|, the results of |Q cannot be explained in terms 
of universal crossover functions. 

On the other hand, it is interesting to note that J~g(t ) 
is non-monotonic. In three dimensions the effect is rather 
small (see Fig. |l|). The minimum value of 7 cff (t) is 
7rfF(i)|min ~ 0.994, so that a non-monotonic behavior 
can hardly be seen. This type of behavior had already 
been observed numerically in two dimensions |2(i(| : in 
this case, however, the effect was much larger. The non- 
monotonicity of 7^ can be predicted analytically by cal- 
culating the first correction to the mean-field behavior in 
the low-temperature phase. One can indeed show that 
7 < ^ F (t ) is increasing for \t\ — > 00. For instance, let us con- 
sider the long-range Ising model introduced in Ref. pl| ] 
and defined by the Hamiltonian 



H = — J(xi — Xj)i 



1*3 1 



(36) 



where J(x) = cR m d for |x| < R m and J(x) = otherwise. 
Setting t=(P c - p)/p c and t = tR 2d /^ d \ one finds for 
2 < d < 4 

7eff = l + C\t\^/ 2 + 0(\t\ d - 4 ), (37) 
where the constants C are given by: 

_ r(3-d/2) 

C " 2<*-V/2(d-2)' m 
C- = -2 2 - d ' 2 [1 - 3 • 2 d - 5 (d - 2)] C+. (39) 

in the high- and low-temperature phases respectively. 
Notice that for d — 3, C~ is negative so that 7^(0 
cannot be monotonic. In two dimensions the subleading 
corrections behave as i _1 logt, again with a negative co- 
efficient. In three dimensions, in Eq. (|37|), logarithms 
appear at next next-to-leading order. 

We have also performed a high-statistics Monte Carlo 
simulation of a d = 3 model of self-avoiding walks with 
long-range interactions pjj . We obtain evidence of the 
existence of a universal critical crossover scaling in the 
large- R limit. The resulting universal curves turn out to 
be in agreement with the two-parameter model predic- 
tions p2,E4j, confirming the general arguments. 
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